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1. Introduction 

Strangeness changing leptonic decays of kaons are an important test of tlie standard model. In 
particular, the Kb decays 

K ^ 7i + l + ViJ = e, n (1.1) 

are a sensitive probe of QCD at low energies In addition, these decays provide for the most pre- 
cise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vus [Q, 0, crucial 
for testing the unitarity of the CKM matrix. The dominant source of uncertainty in the extraction of 
Vus, resides in the experimentally determined quantity |/+(0)yj,.v|. The Ka decay rates were mea- 
sured by BNL-E865, KLOE, KTEV, ISTRA-i- and NA48, for a recent review see [|]. The decay of 
a kaon to a pion, charged lepton and a neutrino is described by the matrix element 

{7i%p')\sY^,u\K+ip)) = ^[ip' + p)Mt) + iP-p')^^f-it)l (1-2) 
where f+{t) is the vector form factor and the combination 

is known as the scalar form factor. The matrix element for the charged pion and the neutral kaon is 
related to Eq. ( [L3| ) by isospin symmetry. For the scalar form factor, expansion at ? = 

A(,)=A(0)(l + Ai-^+l,i.^ + ...), ^4, 

defines the slope Ag and the curvature Ag parameters. Analogously we define the expansion for the 
vector form factor. For a precise determination of |K.s |> it is important to improve the accuracy of 
the parameterizations of the form factors using additional theoretical and experimental information. 
We use inputs from current algebra, perturbative QCD, lattice and chiral perturbation theory to 
provide stringent bounds on the slope and curvature parameters of the parameterizations of the 
form factors. More detailed discussions may be found in 



2. Formalism 

Analyticity is the ideal tool for relating the information from the unitarity cut to the semilep- 
tonic range. The formalism applied in this work (see also |^-^) exploits the fact that a bound on 
an integral involving the modulus squared of the form factors along the unitarity cut is known from 
the dispersion relation satisfied by a certain QCD correlator. For scalar form factor this reads 

. 9x ^ r 7„ n 1 , tlmUo(t) „ , 



.mn„(,)>^^Mf^^|/oMP, (2.2, 
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with t± = {Mk^MjiY- Similar expression, involving the correlator Xi{Q^)' can be written down 
for the vector form factor. We can now use the conformal map f — )• z(f) 

that maps the cut ?-plane onto the unit disc \z\ < 1 in the z. = z{t) plane, with t+ mapped onto z = 1, 
the point at infinity to z = — 1 and the origin to z = 0. This mapping transforms the relations Eqs. 
Q and Q to 

i^^"'^0|g(exp(/0))p</ (2.4) 

where 

g{z)=Mt{z))w{z). (2.5) 

The function w{z) is called outer function and can be calculated analytically in our case. The 
function g{z) is analytic within the unit disc and can be expanded as: 

g{z) = gQ + giz + giZ^ + ■■■ , (2.6) 



and Eq. ( [2. 4] ) implies 

oo 

Igf</- (2.7) 

k=0 

Truncating the above series after finite number of terms gives equations for constraints on the shape 
parameters. Improvement of the bound results if fo{t) is known at a number of real points f/. Using 
a Lagrange multiplier method, we get a determinantal inequality which leads to a quadratic form 
in the shape parameters that is bounded by a known quantity. In order to exploit the knowledge of 
the phase, the dispersion contribution from to ti^ should be subtracted from the pQCD value, 

^ ^"216^1 IWTW ' ^ ^ ^ 



which requires the knowledge of |/o(OP in the region t+ <t < tin- Now Eqs. ( |2!l| ) and ( |2!2| ) can 
be written in the following form 

poo 

/ dtp{t)mt)\^\h{t)\^<l' (2.9) 

for a known weight function p{t), and where the function h{t) is analytic in the f-plane cut for 
t > tin and given by 

h{t) = Mt)m)]-\ (2.10) 
where we implement the phase information by considering the Omnes function: 

/ 1 r d(t') \ 

^(,)=exp(-|.,^j, (2.11, 

where 5{t) is the 7 = 1/2 elastic S-wave Kn scattering phase, in the elastic region and arbitrarily 
Lipschitz continuous above ti^, leading to an extended formalism. This is based on the observation 
first made in ^ which points out that the phase of the Omnes function can compensate for that of 
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the form factor in the region ti^), thereby delaying the onset of the branch point to t-m- The Eq. 



( |2.9| ) can be brought into a canonical form by making the conformal transformation 



2(0 = #-^S' (2.12) 

V ^in ~r V 'in ' 

which maps the complex f-plane cut for t > tin onto the unit disk in the z-plane defined by z = 



Now Eq. (2.S) can be written as 

1 r^^ 

— d0|g(exp(/0))|2</'. (2.13) 
27t Jo 

The function g{z) is now defined by 

g{z)=w{z)(o{z)fom)[0{z)]-\ (2.14) 
where the outer function for the Omnes function given by: 



ft)(z) = exp / at , — -- • (2.15) 

For more information and the review of the formalism, see [0-^]. 

3. Experimental and theoretical information 

We briefly give a description of different inputs used for deriving the improved bounds. 
3.1 QCD correlators 

In the limit >> A^^^^, the correlators Xi (Q^) and Xo{Q^) can be calculated by perturbative 
QCD. Recent calculations to order ^ [l^ give 

1 oc 

1+ = Zi{Q^) = ^^(1 + — -0.062a2 -0.162a,^ -0.176a,4), (3.1) 



lo = XoiQ^) = ^^"l\Z"^ (1 + 1.80a, + 4.65a2 + \5.0a^ + 51 Aa^). (3.2) 



We have omitted the power corrections due to nonzero masses and QCD condensates, as they are 
negligible. We get ;t;i (2GeV) = (343.8±51.6) x lO^^GeV-^ andxo(2GeV) = (253±68) x 10"^. 

3.2 Low-energy theorems 

The symmetries of QCD at low energies are very useful sources of information for our formal- 
ism. The vector form factor becomes equal to the scalar form factor at f = 0. The SU(3) symmetry 
implies /+ (0) = 1 . The deviations from this limit are small due to Ademollo-Gatto theorem. Re- 
cent determinations from the lattice give /+(0) = 0.964(5) [jlTI]. In the case of the scalar form 
factor, current algebra relates the value of the scalar form factor at the Callan-Treiman (CT) point 
Akk = Mj^ — M^. to the ratio Fk/F^ of the decay constants [[l^, p^ ]: 



fQ{AKn) = FK/F^ + AcT. (3.3) 
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In isospin limit, Act 
theory ||l 

to F^/Fk 



3.1 X 10"'' to one loop [Mfl and Act — to two-loops in chiral perturbation 



17| ]. At Ak7c{= —Akk), a soft-kaon result Ql 8| ] relates the value of the scalar form factor 
M-Ak^)=FJFk + Act. (3.4) 



Recent lattice evaluations give Fk/Fji= 1.193 ±0.006 [19, 20 1. A calculation in ChPT to one-loop 
in the isospin limit [14] gives Act = 0.03, but the higher order ChPT corrections are expected to 
be larger in this case. In the present work we use as input the values of the vector and scalar form 
factor at ? = 0. For the scalar form factor we impose also the value fo{AK7c) at the first CT point. As 
discussed in [Q], due to the poor knowledge of Act, the low-energy theorem Eq.(3^) is not useful 
for further constraining the shape of the Kfj, form factors at low energies. On the other hand we 
obtain bounds on Act- 

3.3 Phase and modulus along the elastic region of the cut 

As mentioned in Sec. ^ the bounds can be improved if the phase of the form factor along the 
elastic part of the unitarity cut is known from an independent source. According to the Fermi- 
Watson theorem, the phase of the form factor coincides with the phase of the scattering amplitude 
along the elastic part of the unitarity cut. In our calculations we use below tin the phases from 



[21, 22] for the scalar form factor, and from [23, 24] for the vector form factor. We recall that, 



while the standard dispersion approaches require a choice of the phase above the inelastic threshold 
fin, the present formalism is independent of this ambiguity [||]. Above fin we have taken 5(f) as a 
smooth function approaching n at high energies. The results are independent of the choice of the 
phase for t > fin. We have checked numerically this independence with high precision. 



To estimate the low-energy integral in Eq. ( |2.8| ), we use the Breit-Wigner parameterizations 
of |/+(f)| and |/o(f)| in terms of the resonances given by the Belle Collaboration [ ^ ] for fitting 
the rate of T — >• Knv decay. This leads to the value 31.4 x 10"^ GeV"^ for the vector form factor 
and 60.9 x 10"^ for the scalar form factor. By combining with the values of /+ o defined in Eqs. 



(3.1)-(3.2), we obtain 



/+ = (312 ±69) X lO^^GeV 



T-2 



(192 ±90) X 10 



-6 



(3.5) 



4. Results 

In Fig. [l|, the allowed band for the slope Aq is compared with the experimental determina- 
tions. The slope predicted by NA48 (2007) is not consistent with our predictions. However their 
recent determination confirms our predicted range, see Veltri et.al [Q]. We note that the theoretical 
prediction of ChPT to two loops = (13.9;° j ± 0.4) x 10"^ A^' = (8.0;^;^) x 10""^ is consis- 
tent within errors. For the central value of the slope Aq given above, the range of Aq is (8.24 x 
10""^, 8.42 X 10""*). The same is true for the theoretical prediction Aq = (16.00 ± 1.00) x 10"^ 
Aq = (6.34 ±0.38) X 10"^ obtained from dispersion relations. 

As shown in Fig. |2| for vector form factor, except the results from NA48 and KLOE, which 
have curvatures slightly larger than the allowed values, the experimental data satisfy the constraints. 
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Figure 1: The allowed range for the slope of the scalar form factor, when we include phase, modulus and the 
Callan-Treiman constraint (yellow band). The grey band shows the range without using the Callan-Treiman 
constraint. 
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Figure 2: The best constraints for the slope and curvature of the vector form factor in the slope-curvature 
plane, where the allowed region is the interior of the ellipse. 



We note also that the theoretical predictions A+ = (24.9 ± 1.3) x 10"^ A" = (1.6 ±0.5) x 10"^ 
obtained from ChPT to two loops, and A| = (26.05;°-^}) x 10"^ = (1.29+°;°{) x 10"^ and 
= (25.49 ±0.31) X 10"^ A+ = (1.22 ±0.14) x 10"^ obtained from dispersion relations are 
consistent with the constraint. For precise results, see [Q]. As we mentioned, the same formalism 
can be used to derive regions in the complex plane where the form factors can not vanish. In Fig. || 
we show the region where zeros of the scalar form factors are excluded. If we impose the Callan- 
Treiman constraint, the value of foiAxn), the scalar form factor cannot have simple zeros in the 
range — l.SlGeV^ < fo < 0.93 GeV^. The formalism rules out zeros in the physical region of the 
kaon semileptonic decay. In the case of complex zeros, we have obtained a rather large region 
where they cannot be present. 

For the vector form factors, simple zeros are excluded in the interval — 0.31GeV^ < ?o < 
0.23 GeV^ of the real axis, while Fig. ^ shows the region where complex zeros are excluded. 



6 



Constraining the properties ofK/^ 



Gauhar Abbas 




-5-- 



Ret 

Figure 3: Domain without zeros for the scalar form 
factor: the small domain is obtained without including 
phase and modulus in the elastic region, bigger one 
using phase, modulus and Callan-Treiman constraint. 
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Figure 4: Domain without zeros for the vector form 
factor: the small domain is obtained without including 
phase and modulus in the elastic region, bigger one 
using phase and modulus. 



For more results, see We mention that we do not use as input the soft-kaon theorem, but 
derive bounds on the value at the relevant point. Thus we are able to predict a narrow range 
—0.046 < Act < 0.014 for higher order corrections. 



5. Conclusion 



We have studied the shape of the scalar and vector form factors in the K^i, domain, crucial for 
the determination of the modulus of the CKM matrix element \ Vus\- The results are very stringent in 
the scalar form factor case. The most recent results from NA48 [Q] is consistent with our prediction 
for the slope of scalar form factor and restricts the range of the slope to ~ 0.01 — 0.02. 

Our results show that the zeros are excluded in a rather large domain at low energies, which 
provides confidence in semiphenomenological analyses based on Omnes representations which 
assume that the zeros are absent. Unlike the standard dispersive treatments, our method does not 
require the knowledge of zeros and of the phase above the inelastic threshold. Therefore, model 
dependent assumptions are not necessary. The price to pay is that we only are able to predict 
ranges for shape parameters and zero positions. However, due to the high precision of low energy 
measurements and calculations, the predicted bounds are very stringent. 



7 



Constraining the properties of Kit, 



Gauhar Abbas 



References 

[1] L. -M. Chounet, J. -M. Gaillard and M. K. Gillard, Phys. Rept. 4 (1972) 199. 

[2] M. Antonelli, V. Cirigliano, G. Isidori, F. Mescia, M. Moulson, H. Neufeld, E. Passemar and 
M. Palutan etal, Eur. Phys. J. C 69 (2010) 399 [arXiv: 1005.2323 [hep-ph]]. 

[3] V. Cirigliano, G. Ecker, H. Neufeld, A. Pich and J. Portoles, arXiv: 11 07.6001 [hep-ph]. 

[4] M. Veltri, arXiv: 1101. 5031 [hep-ex]. 

[5] G. Abbas, B. Ananthanarayan, 1. Caprini and I. Sentitemsu Imsong, Phys. Rev. D 82 (2010) 094018 
[arXiv: 1008.0925 [hep-ph]]. 

[6] G. Abbas, B. Ananthanarayan, I. Caprini, I. Sentitemsu Imsong and S. Ramanan, Eur. Phys. J. A 45 
(2010) 389 [arXiv: 1004.4257 [hep-ph]]. 

[7] G. Abbas, B. Ananthanarayan, 1. Caprini, I. Sentitemsu Imsong and S. Ramanan, Eur. Phys. J. A 44 
(2010) 175 [arXiv:09 12.2831 [hep-ph]]. 

[8] I. Caprini, Eur. Phys. J. C 13 (2000) 471 [hep-ph/9907227]. 

[9] R A. Baikov, K. G. Chetyrkin and J. H. Kuhn, Phys. Rev. Lett. 96 (2006) 012003 [hep-ph/0511063]. 

[10] P A. Baikov, K. G. Chetyrkin and J. H. Kuhn, Phys. Rev. Lett. 101 (2008) 012002 [arXiv:0801.1821 
[hep-ph]]. 

[11] P. A. Boyle, A. Juttner, R. D. Kenway, C. T. Sachrajda, S. Sasaki, A. Soni, R. J. Tweedie and 
J. M. Zanotti, Phys. Rev. Lett. 100 (2008) 141601 [arXiv:0710.5136 [hep-lat]]. 

[12] C. G. Callan and S. B. Treiman, Phys. Rev. Lett. 16 (1966) 153. 

[13] R. F. Dashen and M. Weinstein, Phys. Rev. Lett. 22 (1969) 1337. 

[14] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 517. 

[15] A. Kastner and H. Neufeld, Eur. Phys. J. C 57 (2008) 541 [arXiv:0805.2222 [hep-ph]]. 

[16] J. Bijnens and R Talavera, Nucl. Phys. B 669 (2003) 341 [hep-ph/0303103]. 

[17] J. Bijnens and K. Ghorbani, arXiv:0711.0148 [hep-ph]. 

[18] R. Oehme, Phys. Rev. Lett. 16, (1966) 215 . 

[19] L. Lellouch, PoS LATTICE 2008 (2009) 015 [arXiv:0902.4545 [hep-lat]]. 

[20] S. Durr, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, T. Kurth, L. Lellouch and T. Lippert et al., Phys. 
Rev. D 81 (2010) 054507 [arXiv: 1001. 4692 [hep-lat]]. 

[21] R Buettiker, S. Descotes-Genon and B. Moussallam, Eur. Phys. J. C 33 (2004) 409 [hep-ph/0310283]. 

[22] B. El-Bennich, A. Furman, R. Kaminski, L. Lesniak, B. Loiseau and B. Moussallam, Phys. Rev. D 79 
(2009) 094005 [Erratum-ibid. D 83 (2011) 039903] [arXiv:0902.3645 [hep-ph]]. 

[23] B. Moussallam, Eur. Phys. J. C 53 (2008) 401 farXiv:0710.0548 [hep-ph]]. 

[24] V. Bernard, M. Oertel, E. Passemar and J. Stern, Phys. Rev. D 80 (2009) 034034 [arXiv:0903.1654 
[hep-ph]]. 

[25] D. Epifanov et al. [Belle Collaboration], Phys. Lett. B 654 (2007) 65 [arXiv:0706.223 1 [hep-ex]]. 



8 



